Dark-field transmission electron microscopy and the Debye- Waller factor of graphene 
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Graphene's structure bears on both the material's electronic properties and fundamental questions 
about long range order in two-dimensional crystals. We present an analytic calculation of selected 
area electron diffraction from multi-layer graphene and compare it with data from samples prepared 
by chemical vapor deposition and mechanical exfoliation. A single layer scatters only 0.5% of the 
incident electrons, so this kinematical calculation can be considered reliable for five or fewer layers. 
Dark-field transmission electron micrographs of multi-layer graphene illustrate how knowledge of 
the diffraction peak intensities can be applied for rapid mapping of thickness, stacking, and grain 
boundaries. The diffraction peak intensities also depend on the mean-square displacement of atoms 
from their ideal lattice locations, which is parameterized by a Debye- Waller factor. We measure the 
Debye- Waller factor of a suspended monolayer of exfoliated graphene and find a result consistent 
with an estimate based on the Debye model. For laboratory-scale graphene samples, finite size 
effects are sufficient to stabilize the graphene lattice against melting, indicating that ripples in the 
third dimension are not necessary. 

PACS numbers: 61.48. Gh,68.37.Lp,61. 05. J-,63.70.+h,61.72.Dd 



I. INTRODUCTION 

Few-layer graphenes are quasi-two-dimensional, pure 
carbon materials with electronic properties that vary 
markedly depending on the number of layers, how they 
are stacked, and on defects such as grain boundariesPSl 
One of the most effective techniques for determining 
these important structural characteristics is proving to 
be transmission electron microscopy (TEM) PEl For sam- 
ples that can be suspended or mounted on an ultra- 
thin, electron-transparent substrate, TEM has a unique 
combination of advantages. TEM can detect three- 
dimensional ripples and corrugation}^ rapidly map wide 
areas for grain boundaries^lHl and thicknesspl determine 
lattice orientation and mismatchpH21 a nd re solve stacking 
sequence^^ and atomic-scale defect sP 3 -^ However, de- 
spite these capabilities the most effective use of TEM for 
characterizing few-layer graphene has been hindered by 
the lack of a quantitative, analytical model that describes 
how these polydisperse materials scatter electrons. 

We present here a kinematical description of selected 
area electron diffraction from multi-layer graphene, and 
data which illustrate features of the calculation. For most 
other materials a kinematical calculation is not prac- 
tical, since the electron-crystal interaction is so strong 
that multiple-scattering processes cannot be neglected.™ 
These non-linear effects generally require a numerical 
treatment. Simulations based on multislice algorithms 
have been shown to be valuable tools for interpreting 
high-resolution TEM images of graphene, and side-by- 
side comparisons allow thickness and stacking determi- 
nations to be made.^21iiJ However, because graphene 
consists of only a few layers of light atoms, dynamical 



scattering is unimportant in the cases of most interest 
(layer number N < 5). By describing how the diffrac- 
tion peak intensities scale with N, the relatively simple 
kinematical treatment offers direct prescriptions for facile 
grain-boundary mapping and precise layer number deter- 
mination. 

The complete analytic treatment includes a Debye- 
Waller factor, which measures long range crystalline or- 
der in terms of the mean-square displacement of atoms 
from their ideal lattice positions. As was pointed out 
by Peierls and independently by Landau in the 1930's, 
the effects of thermal motion on long range crystalline 
order depend markedly on the dimension of the system 
considered f 12 ^ 13 ! In one dimension long range crystalline 
order is precluded, while in three dimensions there is 
a melting temperature marking the transition between 
phases with and without long range crystalline order. 
The study of phase transitions in the intermediate case, 
two dimensions, has gi ven rise to fruitful ideas concerning 
topological orderfSMlln 1968 Mermin showed that, un- 
der quite general assumptions, crystalline order in two 
dimensions is excluded in the thermodynamic limit 
However, he also noted that the approach to the ther- 
modynamic limit is so slow that it might be irrelevant 
for practically-sized samples. 

In diffraction experiments the Debye- Waller factor de- 
scribes the exponential decay of the peak amplitudes with 
increasing scattering angle. Recent calculations have 
found that graphene's Debye- Waller factor is singular ex- 
cept at zero temperature; 1 - which is to say that the mean- 
square displacements caused by thermal fluctuations be- 
come infinite in the thermodynamic limit. This singu- 
larity is a direct manifestation of the expected lack of 
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long range order for such a two-dimensional system. We 
measure the Debye- Waller factor of a suspended sample 
of mechanically-exfoliated monolayer graphene at room 
temperature, and find it to be more than twice as large as 
the zero temperature expectation. Comparison with an 
estimate of the Debye- Waller factor based on the Debye 
model for the phonon band structure shows that finite 
size effects are sufficient to stabilize the graphene lattice. 



II. THEORY 

The construction of a practical model of image for- 
mation in a modern TEM is a daunting task, since the 
electromagnetic lens parameters are not fixed as in an op- 
tical microscope, and the operative contrast mechanisms 
can vary widely depending on the microscope's imaging 
mode. However, the signals generated in dark-field TEM 
imaging are fundamentally based on electron diffraction, 
and this relatively simple process can be modeled analyt- 
ically. In particular, the selected area diffraction pattern 
generated in a TEM is largely independent of otherwise 
crucial microscope parameters such as spherical aberra- 
tion, electron source type, and objective lens defocus. 
Dark-field images generated with a small objective aper- 
ture have signal intensities that are amenable to a quan- 
titative, analytic description that is almost microscope 
independent. 

We calculate the electron diffraction pattern produced 
by the TEM beam using kinematical scattering theory, 
which has its beginnings in the integral form of the 
Schrodinger equation. The number of electrons dN scat- 
tered into solid angle dft is given by^ 



dN = 
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where / is the beam current, r is the exposure time, A is 
the illuminated area, and e is the electron's charge. To 
find the scattering amplitude /(Ak) we treat the crystal 
potential V(r) as a perturbation and invoke the Born 
approximation. Few-layer graphene represents a nearly 
ideal subject for this approximation, since both carbon's 
small atomic number (Z = 6) and the material's thinness 
act to make the scattered amplitude small. At the end of 
the calculation we will return to discuss the limits of the 
Born approximation's validity. The scattering amplitude 
of a crystal can be writterM 
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Here Ak = k/ — k^ is the difference between the final and 
initial electron wavevectors, and 7 is the Lorentz factor 
(for 80 keV electrons 7 = 1.16). A general position r in 
the crystal has been written r = R; + + r', where R; 
and Rb are the discrete lattice and basis vectors of the 



crystal, respectively. The continuous coordinate r' varies 
over the positions nearer to one atom than to any other. 

We abbreviate the scattering amplitude of Eq. ([2| 
/(Ak) = <S2<S&/ a tom(Afc), which defines the lattice 
sum iS;, the basis sum and the atomic form fac- 
tor / atom (Ak). To evaluate / atom (Ak) we use the 
Yukawa potential to describe the screening of the nu- 
clear Coulomb field by the atomic electrons, taking 
V(r') = (Ze 2 1 1 A-K£or')e~^ r . Here /1 is the inverse screen- 
ing length. Extending the bounds of the integral, which 
properly covers only a single atomic site, to include all of 
space introduces negligible error. In this approximation 
the integral is easily evaluated, giving 
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Scattering in the forward direction is clearly preferred, 
since the atomic form factor ([3| provides a progressive 
suppression as Ak becomes large. 

Previous work on the TEM of few-layer graphene nu- 
merically simulates^ electron diffraction using atomic 
form factors that in turn are numerically estimated^ us- 
ing the relativistic Hartree-Fock atomic wavefunctions of 
a free carbon atom. Expression (3| has the advantage 
of being analytic with a simple physical interpretation. 
For a screening length yT 1 — 33 pm Eq. pi) furthermore 
agrees with the Hartree-Fock calculatiorP^Tfor carbon to 
within 20% for all Ak, and to within 2% relative to the 
Ak = value / at om(0) = 290 pm. 

The lattice sum, or shape factor, Si is a geometric se- 
ries that runs over the N c unit cells illuminated by the 
beam. Performing the sum and squaring gives: 
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where ai and a 2 are graphene's lattice vectors and 
N1N2 = N c = A/A cc \\. The area of a unit cell A cc \\ = 
\/3a 2 /2 with a = |ai| = |a.2 1- In the limit N c — > 00, the 
square of the shape factor is proportional to a comb of 
Dirac (5-functions, since^ 
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For a 1 /mi 2 illuminated area N c ~ 2 X 10 7 , so the infinite 
limit is typically a good approximation here. Writing 
the wavevector difference as Ak = G + g + Ak , where 
G = Vibi + u 2 b 2 is a (discrete) reciprocal lattice vector, 
g = abi + /3b 2 is a (continuous) vector in the Brillouin 
zone surrounding G, and Akj^ is the component of Ak 
perpendicular to the sheet, gives 



|S z (Ak)| 2 ~N c J(a)JC8). 
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Here we have used the orthogonality property a^ • hj = 
2ndij relating the primitive vectors of the direct and re- 
ciprocal lattices. 
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FIG. 1. (Color online) Unit cells and primitive vectors are shown for monolayer, AB-stacked bilayer, and ABC-stacked trilayer 
graphene. On the far right is the graphene reciprocal lattice, with large and small spheres representing strong and weak peaks 
respectively. 



Equation ^ is non-zero only for Ak = G + Ak^, 
so we can restrict our analysis of the basis sum Sb to 
this case. For N parallel honeycomb nets with relative 
displacements Rj, the sum over the basis has the form 
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Here we have chosen to have ai • a2 = a 2 cos(7r/3) and 
hi b 2 = b 2 cos(27r/3), as shown in Fig. [l] The reverse 
convention can be obtained with the substitution v 2 — > 
-v 2 . 

The integers v\ and v 2 index the in-plane diffraction 
peaks at normal incidence. We will refer to peaks that 
occur at the same radius from the reciprocal space origin 
as n th order, where n indexes the number of different 
peak radii occurring between the origin and the peaks 
under consideration. Thus the peaks at |Ak| = || = 
G m in at normal incidence are "1 st order" and those at 
|Ak| = V3G min are "2 nd order". Here S = a/V3 ~ 
142 pm is the carbon-carbon bond length. (See Fig. |4]^a) 
for a graphical representation.) 

For monolayer graphene the sum in Q gives unity and 
the calculation of Sb (Ak) is complete. For multilayer 
graphene the three most symmetric stackings (shown 
in Fig. [I]) are designated AA (simple hexagonal), AB 
(Bernal), and ABC (rhombohedral) , with relative in- 
plane displacements of 0, (j mod 2) x (ai + a2)/3, and 
(j mod 3) x (a! + a 2 )/3 respectively. Layer j + 1 is dis- 
placed by j c ~ 335 j pm in the plane-normal direction 
for j € [0,N — 1]. For the most common case of AB 
stacking, the square of Eq. Q is then: 
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for an even number of layers N. The cases of an odd num- 
ber of AB layers, AA stacking, or ABC stacking present 
no special difficulties and give similar expressions. 

Concluding the second line of Eq. pi is a ratio analo- 
gous to the terms appearing in the shape factor Q, but 



in the case of interest here the number of layers N is of 
order unity so the delta function approximation is not 
appropriate. Equation Q describes how the diffraction 
peak intensities modulate as Ak± is varied, which can 
be accomplished by tilting the sample in the TEM. This 
modulation as a function of tilt angle has been observed 
previouslyP^ and provides one method for determining 
the number of layers present in multi-layer graphene. 

The magnitude |AkjJ for a given diffraction peak 
varies like the sine of the tilt angle, and the sine of the an- 
gle between the tilt axis and the relevant reciprocal lattice 
vector G. For the case of bilayer graphene, Eq. ^ pre- 
dicts that the 1 st order diffraction peaks will show inten- 
sity maxima at angles as small as arcsin(<5/2c) ~ 12° and 
intensity minima at angles as small as arcsin((5/4c) ~ 6°. 
A 2 nd order peak has a local maximum at normal inci- 
dence and can go through a minimum at a tilt angle of 
11° if the corresponding G is orthogonal to the tilt axis. 
All of these features are evident in the data of Refs. HHHl 
Because the 2 nd order peaks are at a local maximum of 
Eq. |8| and the 1 st order peaks are not, the 1 st order 
peaks are sensitive to small variations away from normal 
incidence. This sensitivity to tilt angle can be exploited 
to reveal twinning transitions between AB stacking and 
its mirror image ('AC stacking) in Bernal stacked bilayer 
graphene™ 

Neglecting the curvature of the Ewald sphere, |Akj_| ~ 
at normal incidence and the expression ^ simplifies 
considerably. Table [I] summarizes the possible values of 
|<Sfc(G)| 2 for the three common stacking types. Note that 
for the case of AB stacking the |i5>{,| 2 are identical for sin- 
gle and bilayer graphene in the weak diffraction orders. 
(We enumerate orders of equivalent diffraction peaks ac- 
cording to their distance from the Ak = origin. The 
th , 2 nd , and 5 th orders are strong, while the 1 st , 3 rd , 
and 4 th orders are weak. See Table [i] and Fig. |IJ) This 
degeneracy can have experimental consequences. 

As demonstrated previouslyPS grain boundaries in 
polycrystallinc graphene can be mapped in dark-field 
TEM by placing an aperture in the back focal plane to 
select a certain diffraction peak. If the 1 st order peaks 
are chosen to perform the mapping as in Refs. [5] and |31 
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stacking 


Condition 


weak order 


strong order 


AA 


n.a. 


N 2 




AB 


N is odd 
N is even 


(N 2 +3)/4 
N 2 /4 


4iV 2 


ABC 


N mod 3/0 
N mod 3 = 


1 






TABLE I. Summary of \St\ , the square of the relative struc- 
ture factor, for various graphene stackings. The designations 
'weak' (e.g. 1 st ) and 'strong' (e.g. 2 nd ) order are shorthand 
for (vi + V2) mod 3^0 and = respectively. 



then single and bilayer graphene (and two-thirds of the 
ABC-stacked graphenes) give indistinguishable intensi- 
ties. However, in 2 nd order the simple N 2 scaling breaks 
these degeneracies, allowing the facile visualization of 
changes in the thickness of few-layer graphene. Further- 
more, the larger prefactor in 2 nd order means that these 
peaks also provide more signal in all cases except mono- 
layer graphene. We will see that for monolayer graphene 
the 2 nd order peaks are expected to be 5% dimmer than 
the 1 st order peaks due to the suppression caused by 
the atomic form and Debye- Waller factors, but, because 
of the decreased background at larger Ak, the signal to 
background ratio is still improved in 2 nd order in all but 
the cleanest samples. 

To complete the calculation we integrate Eq. over 
the Brillouin zone surrounding each reciprocal lattice vec- 
tor G, which gives the number of electrons scattered 
into each diffraction peak. It is convenient to do this 
integral in terms of the reciprocal space coordinate pair 
(a, (3), since the shape factor ^ is written in terms of 
(^-functions over these variables. Jacobian transforma- 
tion of the differential solid angle element gives d£l = 
■^(X/ a) 2 dad [3, where A = 27r/fc is the electron wave- 
length and we have taken the variation of the scattering 
angle over a particular Brillouin zone ~ G m i n /(y/3k) to 
be small. For an 80 kV accelerating potential A ~ 4.2 pm, 
so this ratio is ~ 1% or less. Collecting our results for the 
normal incidence case gives, for the number of electrons 
■Wpeak scattered into a peak at G, 
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Here we have used |G| = Ifyvf + v 2 ~ V1V2, defined 
the Bohr radius as = h/ (a c mc) ~ 53 pm, and included 
a Debye- Waller factor e~ 2W . The result ^ has one 
great advantage that might not have been expected at 
the beginning of the calculation, which is that -/V pea k 
only depends on the microscope settings through /, 7, 
and A. Together these variables do no more than modify 
the scale factor multiplying the entire expression. Thus 
Eq. Q does not require special calibration measurements 



and can easily be applied to fit diffraction data acquired 
in any TEM, as we will show later. 

The kinematical diffraction theory developed here is 
not generally considered reliable for crystals^ since the 
cumulative effect of many scattering centers (see the N c 
in Eq. [6]) can make the resultant wavefunction quite dif- 
ferent from the unperturbed wavefunction. For iV-layer 
graphene, however, we can systematically evaluate the 
limits of the kinematical treatment. Equation Q, when 
summed over all v± and t>2, gives the total number of 
electrons scattered from the TEM beam by multilayer 
graphene. The first order Born approximation is valid 
to the extent that this number is small compared to the 
number of incident electrons It / e. The sum over diffrac- 
tion peaks converges rapidly and can be performed nu- 
merically. Taking the Debye- Waller factor e~ 2W = 1 so 
as to include electrons scattered coherently into the ther- 
mal diffuse background,^ we find that a single graphene 
layer scatters 0.45% of the incident beam for an acceler- 
ating voltage of 80 kV. This fraction includes those 'scat- 
tered' into the th order v\ — v 2 = peak, which rep- 
resent almost half (46%) of the total number scattered. 
The total fraction scattered increases roughly like N 2 , as 
shown in Table [Tj For AB stacking the fraction scattered 
is 9% for N = 5 layers and 36% for N = 10 layers. For 
N = 17 layers the first order Born approximation has 
broken down entirely, since the scattered beam is im- 
plied to contain more electrons than the incident beam. 
Thus we expect this kinematical theory to give excellent- 
to-good quantitative agreement for N = 1-5 layers, and 
reasonable qualitative predictions for TV < 10 layers. 

In Eq. (pi) the Debye- Waller factor e~ 2W arises as a 
result of disorder in the lattice, and decreases the number 
of electrons diffracted into a given peak. In graphene 
this factor touches on famous old questions ab out t he 
fundamental stability of two-dimensional crystals In 
an infinite, two-dimensional crystal at finite temperature, 
thermal fluctuations create a divergence in 2W. For an 
anisotropic, layered material like graphene the Debye- 
Waller exponent can be written, 



2W = Akpu 2 



Ak : 



_L«_L> 



(10) 



where the it's refer to the mean-square displacements of 
unit cells from their ideal lattice positions, and the sub- 
scripts p and _L designate the in-plane and normal com- 
ponents respectively. For graphene specifically the mean- 
square displacements have been recently calculated,^ 
with the results u 2 — 16 pm 2 and u]_ = 40 pm 2 at T = 0. 
As a particular example of the classic two-dimensional 
problem, graphene's Debye- Waller exponent is found to 
be singular at non-zero temperatures. A qualitative anal- 
ysis valid for temperatures < IK indicates that finite size 
effects remove the singularity, but no numerical estimates 
more precise than the T = values for the mean-square 
displacements are given 

To provide an analytic expression for the Debye- Waller 
factor that is valid for T 7^ 0, we calculate u 2 within the 
Debye approximation E — hv s k for the phonon band 
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structure. Since the Debye model is best in the infrared 
limit where the singularity occurs, we expect this approx- 
imation to successfully capture the essential physics. For 
Ak ~ G, which corresponds to the normal incidence case 
of most interest here, wc find 



2W = 2 
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Here fe^T is the thermal energy, M is the mass of a 
carbon atom, v„ is graphene's in-plane speed of sound, 
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8^/3ir/95 2 is the square of the Debye wavevector, 
and Xd = fiVsko/ksT is the ratio of the Debye and ther- 
mal energies. A similar ratio x s = hv s k s /kBT is defined 
in terms of the smallest wavevector k s that can be sup- 
ported by the crystal. In the limit that k s — > 0, Eq. ( 11 ) 



illustrates the essential features of the singularity in 2W 
that occurs in two dimensions: u 2 diverges logarithmi- 
cally for T^O, but is finite at T = 0. 

In a real crystal the wavevector k s = 2n/L is lim- 
ited by the size L of the crystal. Since the wavevector 
cannot become arbitrarily small, the divergence in the 
Debye- Waller exponent is effectively regulated. At room 
temperature % ~ 9 and 10 4 for a L — 10 fan 

crystal. Dropping the negligible terms in Eq. (Ill leaves 
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Eq. (12 1 is size and tempcraturc- 
the contribution of zero- 
It is also 



The first term in 
independent, and represents 
point motion to the Debye-Waller factor 
proportional to ft, which highlights its quantum ori- 
gin. Numerically the T = mean-squared displacement 
ul = h/(k D Mv s ) = 16 pm 2 for v s = 2.2 x 10 4 m/s, in 
agreement with the calculation of Ref. [T7J The carbon 
atoms' zero-point motion in graphene is about equal to 
A for 80 kV electrons. 

The second term diverges logarithmically as the crys- 
tal size L — > oo, and increases like TlnT with tempera- 
ture. For experimentally realistic values L = 10 /im and 
T = 300 K, the mean-square displacement u 2 = 44 pm 2 , 
which is to say that the finite temperature correction at 
room temperature is twice the size of the T — value. 
For effective sizes L ranging from 100 nm to 1 cm the 
mean-square displacement ranges from 28 to 69 pm 2 . 

The size of this correction is remarkable: it is not over- 
whelmingly dominant, as might be expected for a for- 
mally divergent term, nor is it small compared to the 
zero point motion, even though room temperature is 
low compared to the in-plane Debye temperature 8c = 
Hv s ku/kB — 2600 K. (In contrast the room tempera- 
ture correction to the T — value of graphene's knock- 
on displacement cross section is tinyP) Furthermore, the 
divergence with crystal size occurs extremely slowly, as 
suggested in Refs. Q£] and [2D While at T = 300 K the 

root mean-square displacement yuj^ is almost 5% of the 

carbon-carbon bond length 5 for a L — 10 /im crystal, it 



is still less than 10% of 6 for an astronomical L = 10 12 m. 
Thus, corrugations or ripples in the third dimensio n 1 5 ^ 6 ! 22 ! 
are not necessary to explain the evident thermodynamic 
stability of laboratory-scale samples of graphene at room 
temperature. Suprisingly, in this case even a mole of car- 
bon atoms, representing a crystal with linear dimension 
L = 126 m, is not enough to approximate the thermody- 
namic limit. 

Equations (11 12 ) imply a melting temperature T m for 



graphene that is weakly size-dependent. While the Lin- 
demann melting criterion can be reformulated in t wo di- 
mensions to circumvent the logarithmic divergenc e 1 23 ^ 24 !, 
such steps are not required for a finite-sized crystal. 
Choosing a standard value^ for the Lindemann param- 
eter C = ^tfJS of 15%, we find T m ~ 3800 K for an 

L = 10 /im graphene crystal. Considering the uncer- 
tainty in £, this reasonable estimate encourages confi- 
dence in the model underlying Eqs. (11 12 ) . 



III. DARK-FIELD IMAGING 

The theory presented has been developed to aid our 
understanding of dark-field TEM images of few-layer 
graphene, and to help refine our procedures for acquiring 
such images. In this section we present TEM images of 
graphene grown by chemical vapor deposition, and im- 
ages of graphene prepared by the mechanical exfoliation 
of natural Kish graphite. Since exfoliated graphene more 
closely approaches an ideal, defect-free system, it pro- 
vides a basis for quantitatively evaluating the successes 
and failures of the theory. CVD graphene samples rep- 
resent a more polydisperse set, and here illustrate how 
dark-field TEM can most effectively characterize vari- 
ables of paramount experimental interest, such as layer 
number and grain boundary structure. 

To prepare CVD samples we first grow graphene on 
copper foil (99.8% pure, 25 /im thick) from methane feed- 
stock, following a procedure described previously.^ The 
growth occurs at 1045° C in a quartz-tube furnace with 
a base pressure of 10 mTorr. To transfer the graphene 
onto TEM grids, we spin-coat a 300 nm thick poly- 
methyl methacrylate (PMMA) support layer onto the 
graphene and etch away the copper foil with a 1.6 M 
solution of FeCl2- After rinsing in de- ionized water the 
graphene/PMMA stack is then scooped onto holey car- 
bon TEM grids and baked at 70° C for 1 hour to remove 
the water. The PMMA is removed in an acetone bath, 
and the grid is finally rinsed with isopropanol. 

We use the method presented previousljEHMIZ] t d e _ 
termine the grain structure of a small suspended region. 
Briefly, we locate the region of interest using bright field 
TEM and acquire the selected area electron diffraction 
(SAED) pattern. Placing an 10 /im objective aperture 
in the back focal plane, we select a certain diffraction 
peak (or peaks) and return the microscope to imaging 
mode. The aperture acts as an electron filter in recipro- 
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FIG. 2. (Color online) The selected area diffraction pattern (a) from CVD graphene suspended over a nominally 3 fim hole 
reveals the presence of two distinct grains. No grain structure is visible in the bright field image (b), where most of the contrast 
is generated by residual PMMA. Dark-field images (c-f) are generated using an objective aperture to separately select both the 
1 st and 2 nd order diffraction peaks, as indicated by dashed and solid circles on the diffraction pattern (a) respectively. All four 
dark-field images (c-f) were acquired with the same exposure time (5 s), and are displayed using identical contrast scales to 
permit fair comparison. In the images (c-d) acquired with the 1 st order peaks the grain boundary is barely visible. The images 
(e-f) acquired with the 2 nd order peaks show much better contrast with obvious moire substructure. Furthermore, image (f) 
clearly reveals a lattice- aligned bilayer island (arrow) that is not visible in image (d). The scale bars represent 500 nm. 



cal space, giving a resultant real-space image with signal 
only from the crystalligraphic orientations corresponding 
to the diffraction peaks selected. 

Figure [2] shows representative data acquired in the 
process of generating dark-field TEM images of CVD 
graphene. The diffraction pattern (a) shows two copies of 
the expected hexagonal structure with a relative rotation 
of 4°, indicating that the field of view contains at least 
one grain boundary. The bright field image (b) shows no 
graphene-specific features, but the dark-field images (c-f) 
map the two grains present and identify each with its spe- 
cific crystallographic orientation. We have found that the 
2 nd order peaks are generally preferred for rapid sample 
characterization. When the 1 st order peaks are selected 
to generate the image the contrast between grains is low, 
to the point that the gross grain structure can barely be 
discerned without image processing. While grain bound- 
aries such as those in (c-d) are perceptible, the lack of 
native contrast makes colorizing or boundary 'gu ides for 
the eye' necessary for many display purposes!^' 

In comparison the grain boundaries seen in images (e- 
f) of Fig. [2j which are acquired using 2 nd order peaks, are 
unmistakable. The analysis of the preceding section ex- 
plains why better contrast is expected from these peaks. 
While the atomic form factor and the Debye- Waller fac- 



tor in 2 nd order are 0.26 and 0.93 of the 1 st order factors 
respectively, the |<5>b| 2 factor is four times larger, for a net 
change of only —5%. However, the diffuse background in 
these images decreases by a factor of ~ 2.8 as the distance 
from the reciprocal space origin increases by a factor of 
v3, leading to the improved signal-to-background in 2 nd 
order. The \Sb\ 2 values given in Table [i] also establish 
that the island visible in Fig. [2Ff ) is AB-stacked bilayer 
graphene, for (assuming N < 4) only AB-stacked bilayer 
has the degeneracy required to make this island invisible 
in 1 st order (Fig. [2](d)) . For samples containing both sin- 
gle and AB-stacked bilayer graphene, comparing the 1 st 
and 2 nd order dark-field images allows faster and more 
conclusive real-space layer mapping than, for instance, 
acquiring a diffraction pattern tilt series for a sequence 
of select areas. 

The 2 nd order peaks also clearly reveal moire substruc- 
ture that is less evident in the 1 st order images (c-d) 
of Fig. |] The upper grain (f) is split by one vertical 
moire stripe about 50 nm wide, and the lower grain (e) 
shows a 3-way intersection of similar stripes, the vertical 
leg of which extends into the moire stripe of the upper 
grain. Previously similar features have been identified 
in Bernal-stacked bilayer graphene, and associated with 
twin boundaries between AB- and AC-stackingPl Follow- 
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ing this suggestion we attribute the stripes seen here to 
narrow regions of bilayer overlap, with the 3- way intersec- 
tion representing the boundary between AB-, AC-, and 
BC-stacked regions. 

CVD graphene shows a wide array of morphologies; a 
given sample can show various numbers of layers, differ- 
ent stackings, and non-trivial grain boundary structures. 
To quantitatively test the analytic diffraction theory de- 
veloped here it is preferable to work with the more homo- 
geneous and defect-free samples that can be obtained by 
mechanically exfoliating natural Kish graphite. We char- 
acterize flakes of mechanically exfoliated graphene using 
Raman spectroscopy, and then transfer these flakes to 
TEM grids for dark-field imaging. The graphene is first 
deposited on 300 nm SiC>2 supported by a silicon sub- 
strate. The flake thickness is determined in a Raman 
microscope (Renishaw In Via) using a 514 nm argon ex- 
citation laser with a ~ 4 /im spot. After characteriza- 
tion, isopropanol is dripped onto the chip until a droplet 
covers its entire surface. A 200 mesh copper TEM grid 
covered by a holey carbon mesh (Protochips C-flat I ^m 
holes, 2 /im pitch) is placed inside the droplet with the 
holey carbon facing the target flake. Once the grid is 
aligned with the target flake the isopropanol is allowed 
to evaporate, which pulls the holey carbon into contact 
with the graphene. Isopropanol is then added to separate 
the grid from the chip, and the grid is allowed to dry in 
air. We find that this procedure reliably provides Raman- 
characterized graphene flakes suspended on a TEM grid. 

Figure [3]Ja) shows a composite image of a graphene 
flake prepared with this procedure. The flake consists of 
regions of varying thicknesses, as is clearly evident from 
intensity variations across the composite image. This 
variation proportional to N 2 is expected, since the dark- 
field image is constructed using the 2 nd order peak in- 
dicated on the diffraction pattern (b). Raman spectra 
(c) collected from regions marked approximately by the 
colored circles in the composite (a) indicate mono-, bi-, 
and tri-layer graphene; the monolayer graphene has an 
approximately 4:1 intensity ratio between the 2D peak 
at ~ 2700 cm -1 and the G peak at ~ 1580 cm -1 , while 
between bi- and tri-layer graphene the ratio goes from 
slightly more to slightly less than one, with characteris- 
tic shoulders on the 2D peak! 221221 A small defect, or '£>', 
peak at 1350 cm -1 is observed because the laser spot is 
wider than the flake. Allowing for the signal from the 
graphene edges, the small size of this peak indicates that 
the sample is high-quality graphene relatively free of de- 
fects. 

One hole in the holey carbon grid, encircled with a 
black ring in Fig. [3^ a), contains suspended mono-, bi-, 
and tri-layer graphene together, along with an actual 
puncture that has no spanning membrane at all. Higher 
magnification images (d-e) of this region, acquired in 
dark-field from the 1 st and 2 nd order spots (designated 
with dashed and solid rings on the diffraction pattern 
(b)), illustrate the relative merits of these two orders for 
rapid sample characterization. In the 1 st order (d) im- 



age the contrast change between the mono- and bi-layer 
graphene is negligible, as expected given the degeneracy 
in the basis sum |6>b| 2 . In the 2 nd order image (e) there 
is a clear contrast progression as the number of layers 
N ranges from to 3. The plot (f) shows normalized, 
average signal intensities collected from regions of dif- 
fering thicknesses (indicated with false color in the in- 
set). A small background, probably due to residue from 
the transfer process, has been subtracted from each in- 
tensity before normalization. The normalized intensities 
collected in 2 nd order quantitatively illustrate the N 2 
scaling predicted for the strong diffraction orders (see 
Table [I]). The dark- field intensities measured using the 
1 st order peaks show the expected degeneracy between 
mono- and bi-layer graphene, along with a signal that 
is approximately three times larger for tri-layer. Again, 
these behaviors are in accordance with the theoretical 
results summarized in Table |U 



IV. MEASUREMENT OF THE 
DEBYE- WALLER FACTOR 

Suspended over the 1 ^m hole marked with a yellow 
triangle in Fig. |3ja) is a single layer of mechanically exfo- 
liated graphene. We perform electron diffraction on this 
region to quantitatively test the analytic model summa- 
rized by Eq. ( 9 ) . A representative diffraction pattern, 
shown in Fig. 4 is acquired with a 0.4 s exposure at a 
camera length of 130 mm using a 2048 pixel x 2048 pixel 
Gatan UltraScan camera with 16 bit digitization. The 
beam blocker prevents measurement of some peaks, but 
roughly 5/6 ths of the total number are visible. The ratio 
of the average intensity of the visible 2 nd order peaks to 
that of the 1 st order peaks is 1.05. Since the expected 
ratios are 0.95, 3.8, 2.9, 3.8, and 3.4 for 1-5 AB-stacked 
layers respectively, this measurement further establishes 
the material as monolayer graphene. 

We now perform a fit to Eq. ^ on 112 peaks extracted 
from the diffraction pattern in Fig. |4|a). Each peak is 
fit individually to a 2D Gaussian with nine free param- 
eters: an amplitude, a 2D linear offset background, two 
center coordinates, two widths, and an angle describing 
the relative rotation between the Gaussian axes and the 
diffraction image axes. All points in the fit region are 
weighted equally A typical width is ~ 2 pixels, the sep- 
aration between peaks is 189 pixels, and a fit region is 
60 pixels x 60 pixels, so the offset gives a good measure 
of the non-Bragg scattering background in the neighbor- 
hood of the specific peak. Near the reciprocal space origin 
the measured peaks have an excess in the wings relative 
to the Gaussian fits, so a peak intensity is calculated by 
summing the counts, minus the offset background, in a 
20 pixel x 20 pixel region about the Gaussian center. 
This procedure gives peak intensities which are insensi- 
tive to the details of the fitting function, except through 
its determination of the background. These intensities 
are plotted in Fig. |4|b) as a function of distance from 




FIG. 3. (Color online) A composite dark-field image (a) of a single graphene crystal on a holey carbon grid is formed using 
electrons from the 2 nd order peak, which is indicated by the solid circle on the diffraction pattern (b). Graphene over a hole, such 
as the one indicated by the yellow triangle in (a), is fully suspended. Raman spectra (c), which were acquired from the regions 
designated by the colored circles in (a) before transfer to the grid, confirm the thickness of each region. Higher-magnification 
dark-field images (d-e) of the region within the black circle in (a) are acquired using the 1 st and 2 nd order peaks respectively. 
The same area is shown again as a false color legend in the plot (f), which depicts the average normalized intensity of each 
disconnected region of (d) and (e) as a function of the square of the number of layers. 



the reciprocal space origin, or, equivalently, the scatter- 
ing angle. The strong peaks have been normalized by 
dividing the integrated intensity by 4. 

We fit peaks as far out as 19 th order, which is at a dis- 
tance 4-\/3G m i n from the origin. The standard deviation 
of the intensities at each diffraction order is taken as an 
estimate of the errors used to provide weights for the fit. 
The fit has three free parameters: a scale factor propor- 
tional to the beam current /, the inverse screening length 
/i, and the mean-square displacement itp, since the other 
parameters in Eq. (I9J) are well known. The final result 
is shown in Fig. |4|b), along with static lattice and zero 
temperature predictions obtained from Eq. (11). 

For the screening length we find = 34 ± 2 pm, 
in agreement with our estimate of 33 pm based on 
the results of the relativistic Hartree-Fock calculation 
in Ref. 19, For the mean-square displacement we find 
— 40 ± 10 pm 2 , where the quoted uncertainty reflects 
a systematic influence of the camera length on the mea- 
sured M 2 that is not presently understood. This value, 
while significantly larger than the T — expectation 
of 16 pm 2 , is consistent with the model described by 
Eqs. (1112) with a characteristic crystal size L ~ 3.5 /Ltm 
at T — 300 K, in accord with the sample dimensions ev- 
ident in Fig. [3]ja). 

Taking the results of the fit at face value, 9% of the 
electrons that would be scattered into diffraction peaks 
by an ideal, motionless monolayer of graphene are instead 



scattered into the thermal diffuse background at room 
temperature, as compared to 5% due to the T = zero- 
point motion alone. While the agreement between the 
observed size of the measured crystallite and the value 
implied by the fit is good, the roles of substrate inter- 
actions, defects and other temperature-independent dis- 
order in the sample, and the actual phonon band struc- 
ture are not yet understood. Future studies that look 
at the Debye- Waller factor as a function of temperature, 
crystallite size, and defect density (e.g. CVD v.s. me- 
chanically exfoliated graphene) will provide further in- 
sight into the nature of long range order in this model 
two-dimensional system. 



V. CONCLUSION 

We have calculated electron diffraction from multi- 
layer graphene using first-principles kinematical scatter- 
ing theory, and presented experimental data from both 
CVD and mechanically exfoliated graphene that illus- 
trate features of the calculation. For a number of layers 
N < 5 the fraction of incident 80 keV electrons scattered 
is less than 10%, indicating that the kinematical theory 
is reliable in this few layer limit. Expressions describing 
the structure factor as a function of N suggest that the 
2 nd order peaks are the most useful for rapid dark-field 
characterization of few-layer samples. For the purposes 
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of mapping grain boundaries, the 2 nd order peaks gener- 
ally give better contrast. For samples containing regions 
of various thicknesses, the 2 nd order peaks give reliable 
thickness contrast oc N 2 and in some cases allow im- 
mediate determination of the absolute number of layers. 
Finally, the complete model includes an analytic calcula- 
tion of the Debye- Waller factor that is valid at T ^ for 
samples of finite size. This part of the calculation exem- 
plifies a famous result for 2D systems that has been taken 
to imply that a 2D crystal is unstable to thermal fluctua- 
tions. Surprisingly we find that, for graphene crystals of 
experimentally relevant sizes at room temperature, the 
thermal corrections to the T = result are several times 
larger than the T = result itself, even though room tem- 
perature is only 10% of graphene's Debye temperature. 
We measure the Debye- Waller factor of a monolayer of 
mechanically exfolliated graphene at room temperature, 
and find it to be more than twice the size of the T = 
expectation. This measurement confirms the importance 
of the finite-size, finite-temperature corrections. We an- 
ticipate that future measurements of the Debye- Waller 
factor in graphene as a function of defect density, crystal 
size, and temperature will clarify the effect of very long- 
wavelength phonons on graphene's strength and stability. 
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